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Abstract 

The lateral Casimir force, which arises between aligned sinusoidally corrugated surfaces of a 
sphere and a plate, was measured for the case of a small corrugation period beyond the applicabil- 
ity region of the proximity force approximation. The increased amplitudes of the corrugations on 
both the sphere and the plate allowed observation of an asymmetry of the lateral Casimir force, 
i.e., deviation of its profile from a perfect sine function. The dependences of the lateral force on the 
phase shift between the corrugations on both test bodies were measured at different separations in 
two sets of measurements with different amplitudes of corrugations on the sphere. The maximum 
magnitude of the lateral force as a function of separation was also measured in two successive 
experiments. All measurement data were compared with the theoretical approach using the prox- 
imity force approximation and with the exact theory based on Rayleigh expansions with no fitting 
parameters. In both cases real material properties of the test bodies and nonzero temperature 
were taken into account. The data were found to be in a good agreement with the exact theory 
but deviate significantly from the predictions of the proximity force approximation approach. This 
provides the quantitative confirmation for the observation of diffraction-type effects that are disre- 
garded within the PFA approach. Possible applications of the phenomenon of the lateral Casimir 
force in nanotechnology for the operation of micromachines are discussed. 
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I. INTRODUCTION 



The Casimir effect [II] is presently well known due to the many potential applications in 
both fundamental physics and nanotechnology. The Casimir force is an extension of the van 
der Waals interaction to larger separations between macroscopic bodies where the relativistic 
retardation of the electromagnetic interaction becomes important due to the finite speed of 
light. The Casimir effect originates from the existence of zero-point and thermal oscillations 
in restricted quantization volumes. The results of extensive experimental and theoretical 
studies on the role of the Casimir force in configurations with idealized boundary surfaces 
(infinitely thin, with ideal- metal boundary conditions etc.) and also between surfaces of real 
material bodies with account of roughness, nonzero skin depth and nonzero temperature are 
presented in Ref. 2] (see also earlier books and reviews p-lO]). Interrelationship between 
experiment and theory in the case of real materials was specially reviewed in Ref. 11 ]. 
In condensed matter physics the Casimir force is taken into account in the investigation 
of various properties of thin films, surface tension, in atomic force microscopy and critical 
phenomena [12]. When the characteristic sizes of microdevices shrink below a micrometer, 
the role of the Casimir and electric forces become comparable. This opens new opportunities 
for the creation of nanoscale devices actuated by the Casimir force [3, [3] • Investigation 
of the combined action of the electric and Casimir forces in microdevices is vital for the 
resolution of issues associated with their stability 



15 



16]. 



The most universally known normal Casimir force acts in the direction perpendicular to 
the interacting surfaces. However, when the material properties of the interacting bodies 
are anisotropic or they are asymmetrically positioned, a lateral Casimir force may exist 
which acts tangential to the surface. Similar to the normal Casimir force, the lateral force 
originates from the modification of electromagnetic zero-point and thermal oscillations by 
material boundaries. For plates made of anisotropic materials, the lateral Casimir force 

n n fin 

and related torques were theoretically predicted in Refs. jl7|, |18| (see also Refs. |4|, 119]). 
For two parallel ideal metal plates covered with uniaxial sinusoidal corrugations of equal 
periods, the lateral Casimir force was predicted [20|, |2l] and calculated [22|, |23j in the second 
perturbation order with respect to the amplitudes of corrugations. In both cases a sinusoidal 
dependence of the result on the respective characteristic angle was obtained. For sinusoidally 



corrugated plates made of real metal described by the plasma model, the lateral Casimir 
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force was considered in Refs. 



24 



25j in second order perturbation theory. In the case of two 



ideal metal parallel plates covered with periodic uniaxial corrugations of rectangular shape, 
some results for the normal Casimir force were found in Refs. j^ . I^ . It is pertinent to 
note that the use of normal and lateral Casimir forces gives the possibility to actuate both 
normal and lateral translations of corrugated surfaces in micromachines by means of the 
electromagnetic zero-point oscillations. The lateral Casimir force might be used to solve the 



tribological problems plaguing the microdevice industry. Specifically, it was suggested 



28- 



30] to use this effect for frictionless transmission of lateral motion. There are also proposals 
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32] to measure the Casimir torque arising between the anisotropic bodies due to the 



lateral Casimir force. 

Both the normal and lateral Casimir forces are in general nonadditive and possess a 
complicated dependence on the geometrical shape of the boundary surfaces connected with 
diffraction effects. In the configuration of a smooth Al coated sphere above a sinusoidally 
corrugated Al coated plate the nontrivial behavior of the normal Casimir force was exper- 
imentally demonstrated 33|. This was done in an additive regime when the sphere-plate 
separation is much smaller than the period of corrugations 23]. Recently (m] the nontrivial 
behavior of the normal Casimir force caused by the complicated geometry of a boundary 
surface was observed in the configuration of a smooth Au coated sphere above a Si plate 
covered with rectangular corrugations (trenches). In this case the nonadditive regime was 
used where separation distances are of the order of corrugation period. The deviations of 



the measurement data from the additive theory were reported 341 ] . These deviations are, 
however, about 50% less than those expected for ideal metal boundaries in the theory taking 
exact account of the geometrical shape 27( . The remaining discrepancy might be explained 
by the role of a nonzero skin depth which was not accounted for in Ref. {27]. Thus, the 
scattering approach for grating geometries taking the exact account of surface geometry at 
zero temperature and describing metallic properties by means of the simple plasma model 
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was found to be in better agreement with the data of Ref. |34j. 
The first measurement of the lateral Casimir force was performed between the surfaces of 
an Au coated sphere and an Au coated plate both covered with aligned uniaxial corrugations 
of sinusoidal shape [37|, |38| . The corrugation amplitudes were A 1 = 59 nm on the plate and 
A 2 = 8 nm on the sphere. The period of corrugations was equal to 1200 nm, i.e., much larger 
than separation distances between the sphere and plate. This means that the experiments of 
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was done in the application region of the proximity force approximation (PFA) 
11] . The lateral Casimir force with an amplitude of 3.2 x 10~ 13 N at the shortest 
separation of 221 nm was found to sinusoidally oscillate as a function of the phase shift 
between the corrugations. The total experimental error of the force amplitude measurement 
at the closest point was determined to be around 19% at the 95% confidence level Q. The 
experimental results were compared with theory using the PFA which should supposedly 
be applicable at separations much smaller than the corrugation period. In so doing the 
dielectric properties of Au were described using the simple plasma model in fourth order 
perturbation theory with respect to a small parameter related to the ratio of the skin depth to 
the separation distance. The experimental results were found to be in good agreement with 
the theory in the limits of the experimental error. It was also predicted theoretically (but 
not observed experimentally due to the use of small corrugation amplitudes and insufficient 
precision) that the lateral Casimir force is asymmetric, i.e., deviates from purely sinusoidal 
harmonic dependence on the phase shift, through the contribution of higher order terms in 
the perturbation expansion with respect to the corrugation amplitudes. 

In this paper we present complete measurement data of the recently performed experi- 
ment on the lateral Casimir force between a sinusoidally corrugated sphere and plate covered 
with Au layers. The much smaller period of corrugations used in this experiment provides 
opportunities to measure deviations of the lateral Casimir force from the theoretical predic- 
tions obtained on the basis of the PFA. Much deeper corrugations on both the sphere and 
the plate and shorter separations where the measurements were performed make it possible 
to demonstrate the asymmetry of the lateral Casimir force, as predicted in Ref. 38(. The 
measurement data are compared with the exact theory describing the Rayleigh scattering of 
the electromagnetic oscillations on the sinusoidally corrugated boundary surfaces of a sphere 
and a plate with no fitting parameters. Only the configuration of a perfect sphere near a 
plane plate is treated using the PFA. At the separations considered, this leads to only a 
negligibly small error of order 0.1% (2), [llj]. The dielectric properties of Au are described by 



means of the generalized plasma-like model [39|, |40[ which takes into account the interband 
transitions of core electrons. The computational results for the amplitude of the lateral 
Casimir force as a function of separation and for the lateral Casimir force as a function of 
the phase shift at different separations obtained using the exact theory at the laboratory 
temperature (T = 300 K) are found to be in a very good agreement with the measurement 
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data. The same data deviate markedly from the theoretical results computed using the PEA. 
This clearly demonstrates the influence of diffraction effects on the lateral Casimir force. 
Some of the results of this paper related to the first set of our measurements were briefly 



published in Ref . [41| . Here, we present additional phase dependences of the lateral Casimir 
force not included in Ref. [4l| for the same corrugation amplitudes, as were used in 41], and 
also the results of further experiments performed with a much larger corrugation amplitude 
on the sphere. We demonstrate that the experimental data of various measurements for 
different samples are in mutual agreement. We also present additional important evidence 
characterizing the quality of corrugated surfaces, procedures of electrostatic calibrations, and 
details of error analysis. Calculations using the PEA are done by means of the exact Lifshitz 
formula rather than with the help of perturbation theory in corrugation amplitudes, as in 
Ref. 4l|. The role of surface roughness is estimated. The details of exact computations of 



the lateral Casimir force using the Rayleigh scattering theory at nonzero temperature given 
in this paper have never been presented previously. 

The paper is organized as follows. In Sec. II we describe the experimental setup and 
sample preparation. Section III contains the description of the procedures of electrostatic 
calibration using both the normal and lateral electric forces. In Sec. IV the measurement 
data and the error analysis are presented. Section V is devoted to the computation of the 
lateral Casimir force between corrugated surfaces using the PFA. In Sec. VI the lateral 
Casimir force is computed in a more fundamental way using the Rayleigh theory. In both 
Sees. V and VI the theoretical results are compared with the experimental data. Section 
VII contains our conclusions and discussion. 



II. EXPERIMENTAL SETUP AND SAMPLE PREPARATION 

A schematic diagram of the experimental setup is shown in Fig. [TJ The experiments 
under discussion are performed using a standard atomic force microscope (AFM) in oil-free 
vacuum chamber at a pressure below 10 mTorr and at room temperature. The observation 
of the lateral Casimir force is done using two test bodies whose surfaces are covered with 
longitudinal sinusoidal corrugations of the same period. The axes of the corrugations should 



be perfectly parallel. Following 



, one can see that misalignment by 1° of the axes of 



corrugations can lead to the loss of any lateral force. As the first test body, we have used 
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a sinusoidally corrugated grating of size 5x5 mm 2 vertically mounted on the piezotube of 
an AFM (see Fig. [TJ). Here, the corrugations have an average period A = 574. 7 nm (i.e., 



more than two times smaller than in Refs. 



37 



38 



in order to achieve the regime where 



the PFA becomes inappropriate for the calculation of the lateral Casimir force between the 
corrugated surfaces). The amplitude of corrugations on the grating A% = 85.4±0.3 nm was a 
factor of 1.45 larger than in Refs. 37J, |38|. A 320 fim long V-shaped silicon nitride cantilever 
of the AFM was first uniformly coated with 40 nm of Al. This modification of cantilever was 
performed to improve its thermal and electric conductivity and prevent deformation due to 
differential thermal expansion in vacuum. 

A polystyrene sphere of a 200 ±4 /mi diameter was placed at the end of the cantilever with 
conductive Ag epoxy. Next, a freshly cleaved mica sheet of 400 /im length, 200 /im width, 
and a few micrometer thickness was attached to the bottom of the sphere also with Ag epoxy. 
Then a second polystyrene sphere with a nominal diameter 2R = 200 ± 4 /im was attached 
to the bottom free end of the mica sheet (see Fig. [1]). After imprint of corrugations on its 
surface (see below), this last sphere was used as the second test body in the measurement of 
the lateral Casimir force. We emphasize that the silver epoxy is rigid at all attachments. The 
first sphere and mica sheet are needed to isolate the laser reflection spot on the cantilever 
tip from the interaction region between the two corrugated surfaces. As a result, the effect 
of scattered light from the top and sides of the corrugated plate is substantially reduced. 
The resulting system (cantilever, mica sheet and the attached two spheres) was uniformly 
coated with a 10 nm layer of Cr and then with a more than 50 nm layer of Au in a thermal 

n 

evaporator. Uniformity of coating was ensured through rotation of the cantilever [42J. The 
Cr layer was used to improve the adhesion of the Au coating and prevent its peeling during 
the imprint of the corrugations on the sphere. 

The lateral Casimir force arises when the uniaxial corrugations on the two test bodies 



are perfectly aligned and have the same period A 



20 



23 



37. 



To ensure that these 

conditions are satisfied, we have imprinted corrugations on the second sphere using the 



grating as a template. In order to obtain deeper corrugations than in Refs. [371, |38], a 
more sophisticated procedure was implemented. Here, the template grating had a 300 nm 
Au coating applied by the manufacturer on top of sinusoidal corrugations made of hard 
epoxy on a 3-mm-thick Pyrex substrate. Note that hard epoxy does not require the use 



of the oxide of Al layer which was coated on the soft plastic grating in Refs. 37], 



381. The 
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smaller period of corrugations used here required the use of a precise stepper motor and 
piezocontrolled imprinting technique. The imprinting procedure was done as follows. At 
first, the second sphere was brought in contact with the grating using a micromanipulator. 
Then a hard flat stylus (see Fig. [T]) was moved in 10 /zm steps using a stepper motor until 
it comes into contact with the other side of the sphere which is now sandwiched against 
the grating. Next the imprinting was done by applying a voltage to the z piezo to gently 
press the sphere between the grating and the hard flat stylus. To obtain deeper imprints 
while preserving sphericity some additional pressure was put on the sphere using the stepper 
motor. Special care was needed at this step. If the applied pressure from the hard flat stylus 
was too much, the resulting rotation of the second sphere will cause the cantilever to break. 
There is less concern on the preservation of sphericity as polystyrene spheres are elastic. 
Then the voltage on the z-piezo is reversed to remove the pressure on the sphere. With 
the help of the same stepper motor the hard flat stylus was gently removed and the sphere 
was then translated horizontally by around 200 /iin to a different position on the grating. 
The last operation is necessary because the part of the grating used to imprint corrugations 
on the sphere might have changed its local amplitude during the imprinting process. This 
translation was performed preserving the orientation of the cantilever. The same imprinting 
procedure, but with application of larger pressure by means of the stepper motor, was done 
with a similar sphere in the second set of our measurements. After the imprint process was 
done, the system was left in vacuum for more than 2 hours to reach a stable equilibrium. 

The corrugations on the grating and on both spheres were examined using an AFM and 
found to be rather homogeneous. This was done after the completion of the experiment. 
An AFM scan of the Au-coated corrugations of the vertical template grating is shown in 
Fig. [2(a). The value of the amplitude of these corrugations A\, also measured using the AFM, 
is indicated above. The variance of the stochastic roughness on the grating was estimated 
in the following way. The topography of the groove from an AFM scan [see Fig. [2(a)] was 
fitted to a sine function. Note that the top and bottom of corrugations are symmetric as 
expected for a sine. This is illustrated in Fig. 12(b) where a typical section of Fig. [2(a) along 
ay — const plane is shown. Then the difference of heights between the groove topography 
and the sine function was calculated at a sufficient number of intermediate points. The 
variance was found using these differences as the root-mean-square deviation between the 
two curves. A total 30 grooves were used to evaluate the average value of the variance 
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leading to 8\ = 13 nm. In Fig. |3]^a) we show an AFM scan of the corrugations imprinted on 
the sphere used in the first set of our measurements. The amplitude of these corrugations 
over a 30 x 30 pm 2 area on the sphere was measured with an AFM to be A 2 = 13.7 ± 0.4 nm 



(i.e., an increase by a factor of 1.7 as compared with Refs. |37l. |38|). Figure HHb) is a typical 
section of Fig. EJ^a) along a y = const plane. It illustrates a symmetry between the top 
and bottom of the imprinted corrugations. The variance of the stochastic roughness of the 
imprinted grooves 82 was found using the same procedure as for a grating with the result 
62 = 2.7 nm. 

In Fig. H] the AFM scan of a 5 x 5 /im 2 area of imprinted corrugations on the sphere 
used in the first set of measurements is demonstrated. Here, the lighter tone shows higher 
points, as indicated in the additional scale. As is seen in Fig. HJ the imprinting procedure 
did not destroy the sphericity of the imprinted surface. An AFM scan of the sphere surface 
with imprinted corrugations used in the second set of our measurements was also performed. 
The amplitude of these corrugations was measured to be A2 = 25.5 ± 0.6 nm (increase by 
a factor of 3.2 as compared with Refs. 37J, |38[). The respective variance of the stochastic 
roughness in this case was 62 = 8.8 nm. Note that the relatively large values of the variance 
of the stochastic roughness presented above (in comparison with the case of flat surfaces 



431 ]) are explained by the fact that the r.m.s. deviation is calculated along the 2-axis, i.e., 



not perpendicular to the corrugated surface at most of the intermediate points. In Fig. [5] 
we show a 5 x 5 /im 2 area of the corrugations imprinted on the sphere used in the second 
set of our measurements. As can be seen in this figure, the grooves are deeper than for the 
sphere used in the first set of measurements, and the sphericity of the surface is preserved. 
The diameters of the spheres with the imprinted corrugations used in the first and second 
set of measurements were measured to be 2R = 194.0 ± 0.3 /j,m and 2R = 194.8 ± 0.3 /mi, 
respectively, using a scanning electron microscope (SEM). The SEM was calibrated with 
a NIST-traceable calibration grating which was independently cross calibrated with our 
interferometrically calibrated AFM. 

The distinctive feature of our setup is that the lateral Casimir force acting along the 
rc-axis tangential to the corrugated sphere and a grating leads to the vertical bending of the 
cantilever. This bending is measured using bicell photodiodes A and B in Fig. [TJ Whereas 
a force acting normal to the test bodies (the normal Casimir force) leads to a torsional 
deflection of the cantilever. The torsional spring constant of the used cantilever K tor was 
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found to be 46 times larger than the bending spring constant ifben (see the next section). 
Because of this the normal Casimir force could lead to only a negligible change in the position 
of the second sphere and in the phase of the corrugations. Thus the setup used is much 
more sensitive to detecting the lateral Casimir force, while simultaneously suppressing the 
effect of the normal Casimir force. 



III. ELECTROSTATIC CALIBRATIONS 

The calibration of the cantilever (i.e., the determination of its spring constants), the 
measurement of the residual potential difference Vq between the sphere and the plate, and 
the determination of the separation on contact zq were done by using both the normal and 
lateral electric forces. We begin with the determination of K tor and Vq by means of the 
normal electric force. For this purpose we have measured the torsional cantilever deflection 
due to the voltages V applied to the grating while the sphere remained grounded. This was 
performed through the measurement of the difference signal between the bicell photodiodes 
shown in Fig. [TJ To determine the parameters of our interest, the measured deflection signal 
at some separation a between the zero levels of corrugations on the grating and on the 
sphere, where the Casimir force is negligible, should be fitted to the theoretical expression 
for the normal electrostatic force between the corrugated surfaces. To attain these ends we 
begin with a brief derivation of such an expression in the sphere-plate configuration which 
is not available in the literature. 

In the configuration of an ideal metal sphere above an ideal metal plate with smooth 



boundary surfaces the exact expression for the electric force is well known 44j . In the range 



of separations from a = 100 nm to a = 6 fim it can be represented 45[ by the following 
polynomial with a relative error less than 10 _4 % 

F^(a) = -27ce (V-V ) 2 <S>(a), 

i= — 1 

where eo is the permittivity of free space and the numerical coefficients q are given by 

c_ 1 = 0.5, c = -1.18260, ci = 22.2375, c 2 = -571.366, 

c 3 = 9592.45, c 4 = -90200.5, c 5 = 383084., c 6 = -300357. (2) 
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The approximate expression for the electric force in the case when the surfaces of a sphere 
and a plate are covered with sinusoidal corrugations can be obtained using the PFA, as 
applied to electric forces. For this purpose we consider the separation along the z-axis 
between any two points z s on the corrugated sphere and z p on the corrugated plate 



z s — z p = a + A 2 sin(27rx/A + ip) — A\ sin(27rx/A). 



(3) 



Here a, as defined above, is the closest separation between a perfectly shaped sphere and a 
plate, and ip = 27txq/A is the phase shift between the corrugations of both bodies. Then, 
according to the PFA, the normal electric force between corrugated a sphere and a plate is 
obtained from Eq. ([1]) by averaging over a period as 



F^(a,p) = -27re (V ~V ) 2 j J dx$( 



Zq Z- 



pi- 



The integral in (Tj0) can be calculated by using the following substitution of ([3]) 



(4) 



z s — z p = a [1 + cos(27rx/A — a)} , 



(5) 



where 



= 0{a, <p) = -{A\ + A\ - 2A X A 2 cos <p) 
a 

tan a = (A 2 cos p> — A\)j [A\ sin <p) . 



1/2 



(6) 



Substituting Eq. (E]) into Eq. 



i^ r (o^) = -27reo(V-y Q 



and performing integration, one obtains 

1 



R 

2a 



a 3 (2 + 3/3 2 ) 



a a 2 (2 + /3 2 ) 
+ c + Cl - + c 2 2R2 +c 3 - 2R3 

a 4 (8 + 24/3 2 + 3/3 4 ) a 5 (8 + 40/3 2 + 15/3 4 ) 
+ C4 — — : h C5- 



+C 6 



a 6 (16 + 120/3 2 + 90/3 4 + 5/3 6 ) 
16^ 



(7) 



It is well known that for corrugated surfaces the use of the PFA may not lead to pre- 



cise expressions for the Casimir force 
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3M25I] . For a static electric force, however, the PFA 
should work much better because in this case the diffraction effects are absent. To verify 
the validity of Eq. (171), we applied a 3-dimensional finite element analysis (FEA) using Com- 
sol Multiphysics |46( to numerically solve the Poisson equation with appropriate boundary 
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conditions in the configuration of corrugated sphere and plate. In the FEA computation, it 
is difficult to use the whole corrugated sphere and plate due to computer limitations. This 
necessitates the use of truncated spheres and parts of the plates. Such a replacement is 
justified as the primary contribution to the force comes from sphere-plate regions which are 
in close proximity. However, the error introduced by the truncation has to be independently 
confirmed to be negligible by varying the size of the regions of the sphere and plate used in 
the computations. 

Computations were performed as a function of the common size / of a truncated section 
of the corrugated sphere and a square corrugated plate. The separation distance and phase 
difference between the corrugations on both surfaces were set to a = 250 nm and tp = 7r/2, 
respectively. The system was enclosed in a grounded rectangular box which represents the 
boundary conditions at infinity. Next, boundary conditions (applied voltages) identical to 
the experimental parameters were assigned to all objects. The rectangular enclosure and 
the sphere were both set to be grounded. The rectangular box was automatically set to 
be at infinity by Comsol Multiphysics. Next the Poisson equation was solved for the given 
boundary conditions. Then the normal electrostatic force was calculated by integrating the 
zz-component of the Maxwell stress tensor over the surface of corrugated sphere and com- 
pared with Eq. (J7J). To make sure that the solution converged, the size of the corrugated 
objects and the number of the surface mesh elements were both varied and the force recal- 
culated. For each I, the number of the mesh elements was increased from one million to 
22 millions until the calculated force converged. The size of the truncated section of the 
corrugated plate and sphere was increased from I = 14 /im to I = 50 /im. The convergence of 
the force was observed for I > 45 /im. The corresponding difference of Eq. ([7]) from the FEA 
results varied from -36% to 2.8% when I increased from 14 jum to 50 /im, respectively, with 
convergence at the largest values of /. Thus, we confirmed from simulations that Eq. (J7|) is 
valid with an error less than 2.8%. It turns out that this error is almost phase independent 
and, thus, does not propagate to the lateral electric force. 

Now we return to the experimental electrostatic calibrations using the normal electric 
force. The deflection signal S^ OI was measured for eight different voltages between -0.52 V 
and 0.47 V applied to the grating at the constant grating-sphere separation a = 1 /im where 
the Casimir force is negligibly small (see below for the determination of absolute separations 
using the lateral electrostatic force). The respective experimental normal electrostatic force 
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is given by 

F£ g (a,<p) = k tm SS B (a,<p), (8) 

where k tor is the normal force calibration constant measured in the units of force per unit 
deflection signal (note that for the attractive force the deflection signal is negative). This 
constant is connected with the torsional spring constant discussed above as 

k toI = K tov m tor} (9) 

where the deflection coefficient m tor is measured in the units of length per unit deflection 
signal. The obtained experimental data for S^ or and for the respective force F^ T from Eq. (jSJ) 
was fitted to Eq. (EJ). The resulting mean values of Vq and k tor found from the fit are 

V = -39.6 ± 1.6 mV, k tor = 7.00 ± 0.08nN/unit S. (10) 

Measurement of the lateral Casimir force as a function of absolute separation requires 
knowledge of the lateral force calibration constant 

^ben = ^ben m ben) (H) 

where m^ en is the bending deflection coefficient. If the separation on contact z is determined, 



abso 
by 



ute separation between the mean values of the corrugations on both surfaces is given 

U] 

Cl = Zq -\- Zpiezo + 'S'nor'^torj (12) 



where z P [ ezo is the distance moved by the plate owing to the voltage applied to the piezo- 
electric actuator, S nor is the photodiode difference signal due to the force (either electric or 
Casimir) . 

The determination of k^en and zq is achieved by measuring the cantilever deflection signal 
due to the lateral electrostatic force which arises when a voltage is applied to the grating. 
The measurements of this signal were performed at small separations from close to Zq to 
zq + 120 nm. Note that calibrations using the lateral electrostatic force were done after the 
measurements of the deflection signal due to the lateral Casimir force are performed, but 
are reported in this section for the benefit of the reader. First, a voltage of 141.456 mV was 
applied to the grating. The sphere was kept at a distance 3.96 nm from zq. The phase shift 
between corrugations was changed continuously at a frequency 0.103 Hz with the x-piezo to 
a maximum translation of 3.3 /im. 
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It is important that corrugated surfaces move parallel to each other as the phase is 
changed with the x-piezo. The grating template was mounted as parallel as possible to 
the x-axis on the AFM piezo. However, experimentally there is always a finite nonzero 
angle between the grating template and the x-axis. When the phase change is introduced 
by moving the x-piezo, this tilt will lead to changes in the separation distance between 
the corrugated surfaces, which in turn would result in a systematic error of decreasing or 
increasing force with phase. This nonparallelity should be corrected before start of the 
experiments using the lateral electrostatic force. A voltage is applied to the plate and the 
lateral electrostatic force is monitored as a function of the phase over many periods of the 
grating. A correcting voltage adapted from that applied to the x-piezo using an adjustable 
voltage divider was synchronously applied to the z-piezo. The applied voltage to the z-piezo 
was changed till the amplitude of the lateral electrostatic force become independent of the 
phase. As an example, in Fig. [6]we show the deflection signal due to the lateral electrostatic 
force versus the phase difference xo between the corrugated surfaces (a) before and (b) after 
the tilt correction was introduced. This was done approximately 25 nm away from contact 
with a 0.2 V applied voltage. It can be seen that the amplitudes of each peak are more 
identical in Fig. Efb). This procedure experimentally verifies that the separation distance 
between the corrugations remains fixed as the phase was changed with the x-piezo. 

The cantilever deflection signal corresponding to the total force acting at the sep- 
aration a = zq + 3.96 nm (sum of the lateral electrostatic and lateral Casimir forces) was 
recorded at 8192 evenly spaced data points. The sphere was moved further away from the 
grating by 5.40 nm to a separation 9.36 nm from z^, and the measurement was repeated. The 
cantilever deflection signal was measured at many other a at the same voltage and also 
for a second voltage of 101.202 mV applied to the grating. To obtain the deflection signal due 
to the lateral electrostatic force alone, it is necessary to subtract from the obtained signal 
5*°* the deflection signal S^t due to the lateral Casimir force (the signal S^t was measured 
before and the respective measurement data are reported in Sec. IV). To perform this 
subtraction, S^ t was fitted to the following sum of harmonics: 



which takes proper account of the fact that the lateral Casimir force is asymmetric (see 
Sec. IV). As the sphere-grating separations in the measurements of the lateral Casimir force 
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(13) 



fc=i 
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and of the total force are not identical, interpolation was used to determine the values of 
Ak(a) at the separations corresponding to the total lateral force. After determination of the 
coefficients Ak(a) from the fit, the signal S^ t was subtracted from the data for the total 
deflection signal S^. The obtained deflection signal S^ t now corresponds to the lateral 
electrostatic force 

F*(a,<p) = k hen S* t (a,<p). (14) 

The experimental force data in this equation were used to determine the calibration param- 
eter fcben and zq from the fit to the theoretical expression for the lateral electric force. 

To derive the analytic expression for the lateral electric force acting between the sinu- 
soidally corrugated surfaces of a sphere and a plate, we _start with the electrostatic energy 
in sphere-plane configuration with smooth surfaces 47. 14 ^ 

E el ' sp {a) = -2ne (V - V ) 2 R 



, R ~ 
c_i in he 

a 



i=0 



i + l\R 



i+l 



(15) 



Here, the integration constant c is equal to [4 



i, n 23 e 

-In 2 -I 1 . 

2 40 126 ' 



(16) 



where < 9 < 1, and q are defined in Eq. (j2J). 

Using the PFA, we obtain the electrostatic energy between sinusoidally corrugated sur- 
faces of a sphere and a plate in the form 



E?{a,<p) = ±j dxE el ^(z s -z p ), 



(17) 



where z s — z p is given by Eqs. (jSi) and (jSJ). Using the definition of the lateral electric 
force [see the formula for the phase shift below Eq. (j3j)] 



Iat 



(18) 



and calculating all integrals in flTTj), we arrive at the result 



F&M = n 2 e (V-V Q )*R AlA2 



a 2 A 



E/aNH-i , 



i=l 



(19) 
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Here, the parameter = (3(a,ip) is defined in Eq. ([6]) and the coefficients Yj = Yi(j3 2 ) are 
given by 

Y 1 =Y 2 = 1, Y 3 = l + ±(3 2 , y 4 = 1 + 

Y, = l + U 2 + U\ Y 6 = 1 + 10/3 2 + h\ (20) 
o o o 

By fitting the force data in Eq. (JHJ) to Eq. ( 1T9|) . the quantities fcben and 2o were found. 
This was repeated for four different electrostatic force measurements and the average values 
obtained are 

z = 117.3 ±3.0 nm, k hen = 1.27 ± 0.06 nN/unit S. (21) 

Using an independent quad-cell AFM measurement of the torsional movement of such a 
cantilever, the torsional spring constant was found to be K tor = 0.28 ± 0.5 N/m. The value 
of the bending deflection coefficient m^, en = 209 ±3 nm/unit S was obtained in the same way 



as in normal Casimir force measurements 
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49| from the change in position of the contact 



point for application of different voltages to the plate during the normal electrostatic force 
measurement. Then the value of the bending spring constant K^ en = (6.1 ±0.3) x 10~ 3 N/m 
was found from Eqs. ( ITTj) and (121 p . By comparing the two spring constants, we find that 
K tor /K^ en 46 which means that our system is really much more sensitive to the lateral 
rather than to the normal Casimir force. 

In recent literature j^c]] there is a discussion that in several measurements of the Casimir 
force the residual potential difference Vq depends on separation, whereas in some other it 
is constant. Because of this, we have performed an additional control measurement of Vo 
at short separation a = 127.3 nm using the parabolic dependence of F{^ t in Eq. f[T9"j) on 
V. For this purpose five different voltages close to the residual potential difference were 
applied to the grating and the cantilever deflection was measured as a function of V leading 
to Vo = — 39.4mV. This is consistent with the result ([TO]) obtained at large separation 
a = 1 (Urn and confirms that in our experiment the residual potential difference is separation 
independent. 

The same calibration procedures, as described above for the first set of measurements, 
were repeated for a second set with a larger corrugation amplitude of 25.5 nm on the sphere. 
Here we obtained a residual potential difference Vq = —28.35 ± 1.15 mV and the following 
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values of the separation on contact and of the lateral force calibration constant 

z = 131.4 ±3.8 nm, k hen = 2.05 ± 0.11 nN/unit S. (22) 

Typically in our second set of measurements with the larger amplitude corrugations the 
calibration errors are a bit larger than in the first, where smaller amplitude corrugations 
were used. 



IV. MEASUREMENT DATA FOR THE LATERAL CASIMIR FORCE AND ER- 
ROR ANALYSIS 



Now we describe how the measurements of the lateral Casimir force were performed. In 
this case the residual voltage Vq, as determined in the electrostatic calibration using the 
normal electric force, was applied to the grating in order to make the electric force equal 
to zero. The x-piezo was used to move the grating along the x-axis and thus change <p. 
The z-piezo, which was independently controlled by an external voltage source, was used 
to change a. The piezo extensions with applied voltage in both directions were calibrated 
using optical interferometry 51] . 

Initially the corrugated sphere was positioned 3.79 nm from the separation on contact 
between the two surfaces zq determined by the corrugations and the highest roughness peaks. 
The thermomechanical drift of the separation distance was measured to be 0.14nm/min from 
the difference in the z-piezo voltage to bring about contact of the two corrugated surfaces 
after a time interval around 30 minutes. A phase shift was introduced by moving the x-piezo 
continuously for a total distance of 3.3 /zm at 0.103 Hz. The photodiode signal corresponding 
to the cantilever deflection was filtered with a low-pass filter with a 30 ms time constant and 
recorded at each of the 8192 points corresponding to x-changes of 0.4 nm. The effect of the 
scattered laser light which would lead to a linear modification of the signal with phase was 
found to be negligible in this experiment. 

Then the separation from zq was increased by 3.6 nm from 3.79 nm to 7.39 nm and the 
deflection signal S^ t was similarly measured as a function of ip and recorded. After this, the 
separation from z Q was increased by 3.96 nm from 7.39 nm to 11.35 nm and the measurements 
repeated. Next S^ t due to the lateral Casimir force as a function of (p was measured at 
separations of 20.05, 32.48, 45.30, 58.01 and 70.86 nm from z . 
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Using the values of kh en and Zo obtained through the electrostatic calibration (see Sec. Ill) 
the measurement data for the deflection signal S^ t were converted into values of the lateral 
Casimir force at every separation a as: 



The resulting lateral Casimir force for a sphere with the corrugation amplitude A 2 = 13.7 nm 
over four corrugation periods is shown in Fig. [7Ja-c) as dots versus the normalized lateral 
displacement x/A between the corrugated surfaces at separations a = 124.7, 128.6, and 
149.8 nm, respectively. Similar results were obtained at all other a listed above. As can 
be seen in Fig. [7J the lateral Casimir force is a periodic function of the phase shift. The 
maximum values of the lateral force, / = max|i^2 t |, decrease with the increase of a. The 
important characteristic feature of the periodic curves in Fig. [7Ja-c) is that they are asym- 
metric, i.e., the dependence of F^ t on x/A is not strictly sinusoidal (note that the sinusoidal 
dependence of the lateral force on the phase shift holds only if the calculation is restricted 
to the lowest order in corrugation amplitudes ~ A\A 2 \ see Sec. V), The true asymmetry 
of the measured lateral force is obvious even without the theory curve. For example, in 
Fig. [7(a) the average shift of each maximum point from the midpoint of two adjacent min- 
ima is (0.12 ± 0.2) A. With the increase of separation, the relative contribution of higher 
perturbation orders in a small parameter AiA 2 /a 2 decreases and the dependence of the lat- 
eral force on the phase shift becomes more close to sinusoidal. The solid lines in Fig. [7J are 
related to the exact theoretical computations. They are discussed in Sec. VI. 

The experimental values of the max|i^2 t | versus separation a are shown in Fig. [8] as 
crosses. The arms of the crosses indicate the total experimental errors determined at a 
95% confidence level using the following procedure. From the data of the phase curves in 
Fig. [7J (and related curves at all the other separations measured) we have calculated the 
mean values / and the variances sj of the quantity / = max \ F^ t \ at each a (see columns 

2 and 3 in Table 1). For the first set of measurements under consideration, the averaging 
was performed over the four periods. The respective Student coefficient is t(i +7 )/ 2 (3) = 3.18 
where 7 = 0.95. Thus, the experimental random error is obtained from the data in column 

3 of Table 1 as 



The main sources of systematic errors in the quantity / are the errors due to the uncertainty 



(23) 



A rand /(«) = s/W3). 



(24) 
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in fcben and due to the resolution of data. According to Eq. (T2TT) . AfcbenAben = 0.047. From 
Eq. ff23l) this leads to the systematic error /(a) = 0.047/(a). The systematic error 
resulting from the resolution of data does not depend on separation A s 2 yst /(a) = 0.035 pN. 
Combining these two errors at a 95% confidence level as random quantities characterized 
by a uniform distribution {2], [32I, 53], we obtain the values of the total systematic error 
given in column 4 of Table 1. Using the rule for the combination of random and systematic 



errors at the same confidence level 12. 52 



531 ]. we obtain the total experimental error A tot f 



as a function of separation presented in the fifth column of Table 1. As can be seen in 
Table 1, all absolute errors are larger at shorter separations and decrease with increasing 
separation. This might be explained by local deviations of the shape of corrugations from 
perfect sinusoidal form which are more influential at shorter separations. The total relative 
error of max|Fj^ t | at a = 121.1 nm is equal to 22.5% (at a 95% confidence level). The 
relative errors at separations a = 124.7, 128.6, and 137.3 nm are equal to 13%, 25.5%, and 
14%, respectively. 

The error in the measurement of the absolute separations, Aa, is a combination of the 
error in the determination of zq, Azq = 3nm, indicated in Eq. ( )2T|) and the error in the 
quantity z piezo + S noT m tor in Eq. f fl2|) . The latter is equal to the half of the first step of the 
z-piezo, A p z m 2 nm. Combining these two errors at a 95% confidence level with the help of 
the same rule, as applied above for the systematic errors {2], [52I 53 K one obtains Aa = 4nm. 
This is a factor of eight improvement as compared with Ref. 38]. Note that the bands 
between the dashed and the solid lines in Fig. M are related to theoretical predictions using 
the PFA and the exact theory, respectively. They are discussed in Sees. V and VI. 

The same measurement procedure, as described above, was applied in the second set 
of measurements for a sphere covered with deeper corrugations of an amplitude A2. The 
resulting lateral Casimir force over five or six corrugation periods is shown in Fig. [9]^a-c) as 
dots versus the lateral displacement at separations a = 134, 156.5 and 179 nm, respectively. 
Similar results were obtained also at separations a = 145.2, 201.6, 224, and 246.6 nm. In 
the same way, as in Fig. [7J the lateral Casimir force at each a is a periodic function which 
is asymmetric (nonsinusoidal). The asymmetry is more pronounced at shorter separations. 
The experimental values of the max \F^ t \ versus separation are plotted in Fig. [10] as crosses. 
Similar to Fig. [8] the arms of the crosses show the total experimental errors of force and 
separation measurements determined at a 95% confidence level. Equation f[2"2"j) and the 
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same statistical procedure, as in the first set of measurements, were used to determine these 
errors with the only difference that for the second set the averaging was performed at all 
separations over five periods [the respective Student coefficient is t(i +7 )/ 2 (4) = 2.78 with 
7 = 0.95]. The total error in the measurement of separation, Aa = 4.7 nm, turned out to 
be a bit larger than in the first set of measurements. The total error of the lateral force 
measurements demonstrates similar irregular behavior on separation distance. This can also 
be explained by the influence of local deviations of groove shape from sinusoidal form keeping 
in mind that with larger corrugation amplitude the role of such deviations should be larger. 
In Table 2 we present the mean values, the variances of the mean, the systematic errors, and 
the total experimental errors (at a 95% confidence level) of the maximum magnitudes of the 
lateral Casimir force at different separations for the second set of measurements. The total 
relative error of the lateral force measurements at separations a = 134, 145.2, 156.5, and 
179 nm varies as 23%, 16%, 14%, and 13%, respectively, all calculated at a 95% confidence 
level. 



V. COMPUTATION OF THE LATERAL CASIMIR FORCE USING THE PROX- 
IMITY FORCE APPROXIMATION 

In this section and in Sec. VI we compare the obtained experimental results for the lateral 
Casimir force with two theoretical approaches applicable in the case of corrugated surfaces: 
the PEA and the exact scattering approach, respectively, with no fitting parameters. As was 
mentioned in the Introduction, the PFA approach was also used to compare with th eory 
the measurement data of the first observation of the lateral Casimir force in Refs. {37! 



However, in those papers the real properties of Au were described using the simple plasma 
model and fourth-order perturbation theory with respect to the relative skin depth at zero 
temperature. Here, we develop a more complete description of the experimental data in the 
framework of the PFA approach based on the Lifshitz formula at the laboratory temperature. 
This will help us to separate diffraction type contributions to the lateral Casimir force which 
are beyond the PFA. 

We start with the Lifshitz formula for the normal Casimir force acting between a sphere 
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and a plate made of a real material (Au) but bounded with perfectly smooth surfaces {2. 8.lll| 

OO nQQ 

F c > sp (a) = k B TRj2 / k ± dk x {In [l-r 2 TM (iZ h k ± )e- 2 * a ] 
1=0 Jo 

+ In [1 - r 2 TE (i^, k ± ) e- 2 * a ] } . (25) 

Here, kg is the Boltzmann constant, £j = 2nkBTl/h with / = 0, 1, 2, . . . are the Matsubara 
frequencies, prime near the summation sign adds a multiple 1/2 to the term with I = 0. The 
reflection coefficients for the transverse magnetic and transverse electric polarizations of the 
electromagnetic field are defined as 

ft 1 ^ £ iQi~ki ct 1 \ rna\ 

t-tmKj, k±) = — , r TE «z, k ±) = — — , (26) 

mi + h qi + k 



lt = ki + ^, kf = kl + e^ (27) 



where 

.2 1.2 

c z " ^ c z 

and = e(i£i) is the dielectric permittivity for the material of the sphere and plate calculated 
along the imaginary Matsubara frequencies. Equation ( )25|) is obtained from the standard 
Lifshitz expression for the free energy in the configuration of two parallel plates by the 
multiplication by 2ttR in accordance with the PFA. Keeping in mind the values of the 
experimental parameters (R ~ 100 /im, a ~ 100 nm), one can conclude that the error 
introduced by the use of the PFA in this case (of about a/R = 0.1% is negligibly 

small. 

Now we consider the sphere and the plate covered with sinusoidal corrugations described 
in Sees. II and III. Within the PFA approach, the approximate expression for the normal 
Casimir force acting between corrugated surfaces of a sphere and a plate can be obtained 
from Eq. (|25|) by replacing a with z s — z p defined in Eq. ([3]) and averaging over the period 
of corrugations A. By expanding also the logarithms in Eq. ( )25l) into a series, we arrive at 

OO 



^ r . , k B TR^ 1 f°° , „ 

— l—n •'O 



n=l 1=0 

x [r&(ifc, fcj.) + r&(i& AU.)] e- 2 * na (28) 

-A 

X / ^ x e ~2qin[A 2 sm(2TTx/\+ip)-A 1 sm(2TTx/A)] 







From this expression it is simple to obtain the Casimir energy in the configuration of a 
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sphere above a plate both covered with sinusoidal corrugations 

poo 

E c (a,<p) = / dzF^cp) 

J a 



k R TR^ 1 ^/ f°° ki dk\ 



2A ^ n 2 ^ J Q q t 

n=l Z=0 1/0 

x [r^(i&, + r&(i&, fcj] e- 2 ""- (29) 

r-A 

X / dx e - 2l Jr«[^2sin(27rx/A+v)-Aisin(27rz/A)] 







Introducing the dimensionless variables 

2< 



y = 2ag,, = — (30) 

c 



and performing the integration with respect to x, one can rearrange Eqs. ( 1281) and ( 1291) to 
the form 



=--£5- E-E / 

71=1 Z = 

x Km«^) + r&(iG, v)] e"" y I o (n0y), (31) 



n=i z=n 17 « 



71 = 1 z=o 

x [r* (<0, y) + r* (iO, y)] e""" / (^y). (32) 

Here, Iq(z) is the Bessel function of an imaginary argument and /3 is defined in Eq. ([6]). 

The general expression for the lateral Casimir force can be obtained from Eq. ( )32|) in the 
same way as the lateral electric force in Eq. ( II 8p . After performing the differentiation, we 
rearrange the result to the form 



F*(a,<p)= 2q3A/3 ^2^-1, J (i y d y 

x [r 2 T n M «i, V) + y)] e~ ny h(n(3y). (33) 

This expression is convenient for numerical computations. 

Equations ( 13TI) and ( 1331) are generalizations for the case of real materials of the previously 
known expressions obtained at T = for ideal metal corrugated bodies. In this case r^ M = 
r TE = 1 an d one obtains 

fc„D 00 1 POO 

FZ(a, V ) = -—J2- y'dye-^Wy), (34) 

oTVCL Ti In 

71=1 J{i 

^t(a^) = ^^sin^^i j( tfdye^hWv). 
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Introducing the new variable v = ny and performing the summation, we arrive at 

7i 3 hcR 



(a,tp) 
a,<p) ■ 



720a 3 J() 



aW e- v I (/3v) 



(35) 



^KcRA x A 2 . [°° 2 _ 



After performing the integration in these equations, the final result is 

n 3 hcR 2 + f3 2 



720a 3 (l-/3 2 ) 5 /2' 
7i 4 hcRAiA 2 



(36) 



360a 4 A(l - /3 2 ) 5 / 2 
The last expression for the lateral force was obtained in Refs. 



sin if. 
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38| with corrections due 



to the skin depth. In the first perturbation order in A\A 2 (i.e., for (3 = 0) this expression 



was also used in Refs. 
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25]. 



Equations ( I3TT) and (1331) take into account the sphericity of the second test body and 
sinusoidal corrugations on both bodies in the framework of the PFA and the nonzero skin 
depth in the framework of the Lifshitz theory. To perform numerical computations, one 
should use a model to represent the dielectric permittivity. At the separations considered 
(from 100 to 250 nm) the thermal effects do not play any role in any of the theoretical 
ap pro aches and reliable results can be obtained using the generalized plasma-like model 
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9j 



it ^ uj + a + 7i a 



(37) 



where u p = 9.0 eV for Au is the plasma frequency. The values of the oscillator frequencies 
Uj, relaxation parameters jj and the oscillator strengths Qj leading to approximately the 
same results for Im.e{uj) as those obtained from the optical tabulated data 57J over the 



4Q|. 



frequency range from 2 to 25 eV are listed in Refs. 

An important factor that is not taken into account in Eq. (1331) is the surface roughness. 
Stochastic roughness covering the corrugated surfaces in both the sphere and the plate is 
clearly seen in Figs. El|3] and its measured variances are presented in Sec. II. In view of its 
stochastic (nonperiodic) character, surface roughness cannot influence the phase dependence 
of the Casimir energy in the configuration of corrugated surfaces. However, it contributes 
to the Casimir energy, and, thus, to the lateral Casimir force, through a phase independent 
correction factor which depends only on separation. It has been found 0, Q] that the 
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correction factor to the Casimir force due to stochastic roughness in the configuration of 
a sphere above a plate with no corrugations is given by 



X2 — I — s.2 (p , x2\2 

VF(a) = 1 + 6^±p + 45 ^7 2) . (38) 



Now we take into account that in the framework of the multiplicative approach the correction 
factor due to surface roughness is the same for real and ideal metals. Keeping in mind that 
the Casimir force between a perfectly shaped ideal metal sphere and a plate is F(a) ~ 1/a 3 , 
we obtain the correction to the Casimir energy by integrating the Casimir force with the 
inclusion of roughness, rj F (a)F(a), with respect to a 



A2 _|_ X2 (X2 , fi2\2 

VE(a) = 1 + 3^±p + 15 [dl + / . (39) 



Note that for the experimental variances (see Sec. II) the fourth order term in Eq. (|39|) is 
negligibly small. Thus, for the first set of measurements at the shortest separation a = 
120 nm the second order term contributes 3% of the force, but the fourth order term only 
0.2% of the force. For the second set of measurement with deeper corrugations on the sphere 
the contributions of the second and fourth order terms at separation a = 134 nm are 4% and 
0.3%, respectively. 

To obtain the correction factor to the lateral Casimir force due to the surface roughness 
with account of the sinusoidal corrugations on the sphere and on the plate, we should replace 
a with separation z s — z v between the corrugated surfaces in Eqs. ([3]) and ([5]) and perform 
the averaging with respect to the phase shift and the period 

Vcorr (a) = 1 + 3(5! + $1) ( , \ , (40) 

\ [z s z p ) I 

1 \ 1 f 2w , f A , 1 

dip / dx- 



z p ) 2 I 2vrAa 2 J r J [1 + cos(27ra;/A - a)} 2 ' 

Calculating both integrals in the second equality in Eq. f l4"0"j) and substituting the obtained 
result into the first equality, we arrive at 

5 2 + 5 2 



?7corr(a) = 1+3 



a(a + A 1 -A 2 y/ 2 (a + A 2 -A 1 y/ 



2 



r ,1 1 4A X A 2 . 

1 2' 2' '(a + A 1 -A 2 )(a + A 2 -A 1 ) ' ' 



where F(a, (3; 7; z) is the hypergeometric function. 
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Substituting the values of A* and Si from the first set of measurements (see Sec. II) in 
Eq. ()4ip . we find that at separations a = 120, 150, and 200 nm the correction factor due 
to surface roughness takes the values r] COTT = 1.054, 1.029, and 1.015, respectively. For the 
second set of measurements (i.e., using A 2 and 5 2 in Sec. II) the correction factor in Eq. ( 14 ip 
takes the values fj COTT = 1.060, 1.041, and 1.021 at separations a = 130, 150, and 200 nm, 
respectively. It is seen that even for the second set of our measurements, where the surface 
roughness on the sphere is larger than in the first one, the roughness correction is rather 
moderate and does not exceed 6% of the lateral Casimir force. 

In Fig. |HJ (the lower dashed line) we present the computational results for the maximum 
values of the lateral Casimir force as a function of separation in the framework of the PFA 
[i.e., using Eqs. (I3"3"j) and ( 13"T|) ] for the first set of our measurements. This line does not 
take surface roughness into account. The upper dashed line in Fig. [HJ presents the results 
for max \F^ t \rj COII as a function of separation, i.e., taking into account the correction due to 
surface roughness. Keeping in mind that the multiplicative approach provides the means 
to estimate the effect of roughness rather than to calculate it precisely on the basis of the 
fundamental theory, we consider the band between the dashed lines in Fig. [8] as a theoretical 
prediction given by the PFA approach. As can be seen in Fig. [HJ the experimental data 
indicated as crosses are inconsistent with the prediction of the PFA approach, thus revealing 
the role of diffraction-type effects on the corrugations of relatively small period used in this 
experiment. 

In Fig. Unia similar comparison between the experimental data and theory using the PFA 
approach is done for the second set of our measurements with corrugations of larger ampli- 
tude. The lower dashed line is computed using Eqs. (133]) and (137|) with the surface roughness 
disregarded. The upper dashed line is obtained by the multiplication of the computation 
results by an additional factor fj COTT computed using Eq. (j4Tj) . Again the theoretical band 
between the dashed lines is inconsistent with the experimental data shown as crosses due 
to the neglect of diffract ion- like effects. In the next section we present a more fundamen- 
tal theoretical approach taking the diffraction effects into account which brings theory into 
agreement with the experimental data. 
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VI. COMPUTATION OF THE LATERAL CASIMIR FORCE USING THE 
RAYLEIGH THEORY 



The PFA approach used to calculate the lateral Casimir force in the configuration of a 
sphere above a plate with corrugated boundaries has an uncontrolled error which increases 
with decreasing period of corrugations. As was mentioned in the Introduction, in such cases 
a more fundamental theory is desirable which takes exact account of diffraction effects. 
Here, we present such a theory, based on the scattering approach, for the configuration of 
two plates covered with periodical corrugations of arbitrary shape but common period. This 
theory can be applied to our experimental configuration by using the PFA for the transition 
from the configuration of two plane plates to a smooth sphere above a plate. According to 
our discussion in Sec. V, the error in computations introduced by the use of the PFA for 
this restricted purpose alone, is extremely small and quite satisfactory for the needs of our 
experiment. 



The scattering approach started with the work 58] where it was applied to obtain the 



Lifshitz formula and Refs. 
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60] on the multi-scatte ring ex pansion. On this basis the 
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631 ] and was applied to several 



multipole scattering technique was developed in Refs. 
configurations with curved boundaries for scalar and electromagnetic cases. Alternative 
techniques within the scattering approach to the Casimir effect were presented in Refs. 



671 ] . For objects o 



used in Refs. [54 



sp 
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herica 



and cylindrical shapes the multipole scattering technique was 
68M70| . This technique works well for large separations between the 
objects and shows poor convergence when the separations between the objects are small. At 
small separations the leading asymptotic expansion of the Casimir energy at zero tempera- 



ture within the scattering approach was developed in Ref. [7l|. The technique developed in 



Refs. 
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36] uses another basis, the Rayleigh basis, which is a generalization of the plane 



wave basis. The Rayleigh basis is a natural choice for grating geometries. 

We consider two parallel 3-dimensional longitudinal (along the y-axis) periodic (along 
the x-axis) dielectric (metallic) gratings of arbitrary form separated by a vacuum gap so 
that they form a waveguide and one grating is located above the other (we assume that the 
edge of the corrugation region of each grating is perpendicular to the z axis). The periods 
of both gratings A are equal. In general, the top and bottom gratings can have different 
shapes and different dielectric (metallic) properties. In Fig. [TT]the height of the corrugation 
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region is equal to h for the bottom grating. The meaning of the lateral displacement x$ 
becomes obvious from the comparison of Fig. [11] with Fig. [12] and Fig. [131 We suppose that 
the space between the two gratings is vacuum with e — // = 1 and we assume \i — 1 inside 
the medium. 

The physical problem is time and y invariant, so that the particular solutions for the 
electric and magnetic fields can be written in the form: 



Ei(x,y,z,t) = Ei(x,z)exp(ik y y - iurt), 
Hi(x,y,z,t) = Hi(x, z) exp(ik y y — iut). 

The solutions of Maxwell equations should satisfy the quasi-periodicity conditions: 

Ei(x + A,z) = e ik * A Ei(x,z), 
Hi(x + A,z) = e ikxK HAx,z). 



(42) 



(43) 



In every Casimir problem one needs to determine the complete basis of solutions. Let 
us suppose that the top grating is absent. We consider a generalized conical diffraction 
problem (k y ^ 0) on the lower grating for the incident wave with the x-component of the 
wave vector equal to some fixed value which we denote k x . The longitudinal components of 
the electromagnetic field outside the corrugated region (z > h) may be written by making 
use of the Rayleigh expansion 



Ey(X,Z 



Hy(X,Z 



72] for an incident monochromatic wave: 
) = 1^ exp(ik x x — ifl^z) 

fox 

OO 

+ 2J exp(ia n x + if3^z), 

n=— 00 

) = li h) exp(ik x x-ip {1) z) 

OO 

+ R { njexp(ia n x + i(3 { n l) z). 



(44) 



Here 



with 



oc n = k x + 27m/ A 



2tt 

T 



Ak x 



2tt 



(45) 



(46) 
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where [r] is the integer part of the number r. From this it follows that < k x < 2tt/A. The 
other notation in Eq. ( T4"4"j) is as follows: 

P {1)2 = " 2 ~kl- kl, ^ 2 =^-kl- al (47) 

The quantities I- 6 '^ and in Eq. ()44p are the Rayleigh coefficients. [Note that with 

notations (143]) - (147]) it holds k x = a p , 0^ = /3jp .] The solution ( f44|) is also valid outside 
any periodic structure in x direction. All other field components can be expressed in terms 
of the longitudinal components E y ,H y by using standard formulas as done in waveguide 
theory. This can be done since the factor exp(ik y y) is conserved after the reflection of the 
electromagnetic wave from the grating. 

At z = the solution has to satisfy the expansions 

oo 

E y (x,z) = T$exp{ta n x-tf3i 2 h), (48) 

n=— oo 

oo 

H y (x,z) = ^2 exp(ia n x - i(3 (2) z), 

n=— oo 

where 

/3i 2)2 = a, 2 - kl - al (49) 

which are valid for z < 0. The coefficients Tnp h ^ are called the transmission matrix coeffi- 
cients. 

The E y component of the electromagnetic field in the region < z < h is defined as 
follows: 

oo 

E y (x,z)= E^(z)exp(ia n x). (50) 

n=— oo 

The other components of the electromagnetic field in the region < z < h are defined in anal- 
ogy. Inside the corrugation region < z < h it is convenient to rewrite Maxwell equations in 
the form of the first order differential equations, dA/dy = M(z)A, where M(z) is a square 
matrix of dimension 8iV + 4, A T = (Ey . . . E~ N , E% . . . E~ N . . . H~ N , H* . . . H~ N ) and 
2N + 1 is the number of coefficients considered in the expansion (loT)|) for E y and similar 
expansions for E x , H y and H x . 

For a rectangular grating the matrix M(z) is constant (independent of z). From Maxwell 
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equations 



we get 



OEy 

ik y E x - = -iuH z , (51) 
ik y H x - = iueE z 



E * = TZ ( ik v H * ~ ) > ( 52 ) 



itue \ dx 
H z = -— IKE* 



l - ■ v y J — / x <-. 

tco \ OX 

Substituting into Eq. ( 152]) the expansion ( 150]) and analogous expansions for all other com- 
ponents of E and H, we arrive at 



Ez = E(— ) ('/.„//;." -m, !; //; 

— \IUJ£ / 



7(1 



ra— m 



H n z = -l(ik y E^-ia n E^), (53) 

where (0) n _ m is the Toeplitz matrix. 

The following equations follow from the remaining four Maxwell equations: 

dE™ 



fh - ik y E2-iu;H2, 
dE n 

= ia n E? + iujHl (54) 
^ -- tk y H: + J2(^e) n - m E™, 



dz 
dR 



= ia n H n z - ^(iue) n „ m E™. 



One can substitute Eq. ( 1531) into Eq. ( 1541) and obtain a system of the first order differential 
equations for the Fourier components of the electromagnetic field E y , E™, H y , H™ in the 
region < z < h. 

Now we have to determine the Rayleigh coefficients nip , R^J for the specific periodic 
geometry profile. One can determine these coefficients by matching the solution of equations 
inside the corrugation region < z < h with Rayleigh expansions (1441) at z = h and 
expansions ( 148]) at z = 0. This can be done by imposing the continuity conditions on each 
Fourier component of the fields E y ,E x , H y , H x at z = and z = h. 
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There is no separation of the TE and TM modes for the general case. That is why the 
reflection matrix R^ t for the reflection from the bottom grating can be defined as follows: 

p( e ) /r(e) _ r Ah) _ n \ r>(e) / Ae) _ n Ah) _ c 

i? bot (fc±,w)= | ... _ ... ^ ... ,, , h . _! . (55) 



pO) /r(e) _ r t(7i) _ n \ dO) /r( e ) _ n rC 1 ) _ A 

- rt "3<?3l-'P — °P<?3' 1 P — U ) Jr Hl4,q4\ 1 P — U 1 1 P — U \ 



PQ4 ) 



To obtain the Casimir energy we need to determine the eigenfrequencies of all the normal 
modes of the electromagnetic field between the two periodic gratings. These eigenfrequencies 
can be summed up by making use of an argument principle, which states: 

X f d 

2wi T ^^du ln f^ duj = Yl ~ ( 56 ) 

Here, u are zeroes, are poles of the function f(co) inside the contour of integration, 
and the degenerate eigenvalues are summed over according to their multiplicities. For the 
Casimir energy we have 4>(u) = frw/2. The equation for eigenfrequencies of the corresponding 
problem of classical electrodynamics is f(oo) = 0. 

Consider first the plate-plate geometry when the two dielectric parallel slabs (slab 1: 
z < 0, slab 2: z > L) are separated by a vacuum gap (0 < z < L). In this case the TE and 
TM modes are not coupled. The equation for the TE eigenfrequencies is: 

f{u) = 1 - r« bot (fc ± , u;)4| top (/c ± , w) = 0. (57) 

Here bot (fcj_, uS) is the reflection coefficient of a plane wave moving down which reflects 

(2) 

on the dielectric surface of the slab 1 at z = 0, while r T ^ top (k±,u) is the reflection co- 
efficient of an upward moving plane wave which reflects on the dielectric surface of the 
slab 2 at z = L, and k± = (k x ,k y ). One can deduce from the Maxwell equations that 
r TE,to P ( fc ±> w ) = r TE,bot( k ±' UJ ) e M 2 i k zL) [here, 

r TE,bot(^-L; u ) i s the reflection coefficient 
of a downward moving TE plane wave which reflects from the dielectric slab 2 temporarily 
located at the position of the slab 1, i.e. at z < 0]. From ( |57|) and the analogous equation for 
the TM modes one immediately obtains the Lifshitz formula by making use of the argument 
principle ( |56l) . 

For two periodic dielectrics or gratings separated by a vacuum gap one has to consider 
the reflection of downward and upward moving waves from the unit cell < k x < 2tt/A. 
Imagine that we removed the top grating in the system. Then the reflection matrix of the 
downward moving wave is defined as R^t- Imagine now that we remove the bottom grating 
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(2) 

in the system. Then we denote the reflection matrix of the upward moving wave as Rl L. The 
reflection matrices R^ t , -Rtop depend on the wave vectors of the incident waves, parameters 
of the gratings and the mutual location of the gratings. The equation for normal modes 
states: 

#ba t ( fe ±> ^i)#top( fc ±, u u L, x )ipi = ipi, (58) 

where ipi is an eigenvector describing the normal mode with the frequency U{. Instead of 
Eq. (!57l) one obtains: 



det 



I-RW(k ± ,u)B™(k_ L ,u,L,x ) =0. (59) 



For every k x , k y the solution of fl59l yields possible eigenfrequencies Ui for the solutions 
of Maxwell equations that should be substituted into the definition of the Casimir energy 

(2) 

Suppose that the reflection matrix i?^ ot for the reflection from the fictitious (imaginary) 
grating located as in Fig. [12] is known in the coordinates {x,z). Performing a change of 
coordinates z = —z\ + L, x = x\ — x (x < A) in (PHI) , it is possible to obtain the matrix 

(2) 

Rl op for the reflection of upward waves from a grating with the same profile turned upside- 
down and displaced from the lower grating by Ax = Xq, Az = L (see Fig. [T3|) . It follows 
that 

R®(k ± ,i£,L,x ) = Q*(x )K(k x ,iC,L)R^ t (k ± ,iC)K(k ± ,iC,L)Q(x ), (60) 

where R^ t (k±, i£) is the reflection matrix of the downward moving waves from the grating 
in the system of coordinates (x, z) depicted in Fig. [12j Here, K(k^,il;, L) is the diagonal 
2(2iV + 1) matrix of the form: 



K{k ± ,i(,L)= \, (61) 




with the matrix elements exp[— £ 2 + ky + (k x + 2irm/A) 2 ] (m = —N, . . . , N) on the 
main diagonal of the matrix G±. The lateral translation 2(2A^ + 1) diagonal matrix Q is 
defined as follows: 

(g 2 o \ 

QM = , (62) 

\0 G 2 J 

with matrix elements exp(27rzmxo/A) on the main diagonal of the matrix G 2 , where m = 
—N , . . . , N. 
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The summation over the eigenfrequencies is performed by making use of the argument 
principle ( 156]) . which yields the Casimir energy of two parallel gratings on a unit cell of a 
period A and unit length in the y direction: 



E c (L,x ) 



he A [™ f 2n/A 71 , , 

a£ I dk v I dk r In det 



(2vr) 3 



y 

oo Jo 



I - RV(k ± ,iZ)BW(k ± ,iZ,L,x ) . (63) 



Here, as was defined in Sec. II, (p = 27txq/A and xq is a lateral displacement of the two 
gratings. This is an exact expression valid at zero temperature for two arbitrary parallel 
gratings with equal periods A separated by a vacuum gap. 

The Casimir free energy on a unit surface J 70 for the system of two gratings can be 
written as follows: 

k B T 



7T- 



E 

1=0 



X) pit/ A 

dky / dk x In det 
Jo 



(64) 



This formula is valid for an arbitrary profile and an arbitrary dielectric permittivity of each 
grating. 

When one grating has a curvature of the sphere of the radius R 3> A, the lateral Casimir 
force in this system F^ t can be obtained by combining the use of the PEA for a sphere-plate 
configuration and the exact formula for the free energy of the two gratings fl64l) : 

,dF c {z',<p) 



dip 



(65) 



For the sinusoidally corrugated sphere and plate used in our experiment L = a + A\ + A2 
and the lateral Casimir force is given by 

Ak B TR 



A 



x In det 



so 00 t Q roo rn / A 

dz'y —— / dk„ I dk q 
U^h V Jo 



I - R^ t (k ± , t^)R^ p (k±, z' + A 1 + A 2 , ip) 



? (2) 

Hop 



(66) 



Different methods were developed in grating theory starting with the pioneering work of 



' Jlayleigh 



721 ] . For our case of sinusoidal gratings we used the so-called differential method 



731 ] . The integration of the first-order ordinary differential equations (I54p with the dielectric 



permittivity of the generalized plasma-like model fl3T|) in the corrugation region of each 
grating was based on an explicit Runge-Kutta (4, 5) formula, the Dormand-Prince pair, and 
was performed using the Matlab package. Then the reflection matrices R^ t and R^ p were 
computed. The lateral Casimir force F^ t was computed using Eq. f )66|) for the experimental 



32 



parameters of the first and second set of measurement. Every Matsubara term in Eq. (|66|) 
was evaluated with a precision of 1.5% which is also the computational precision of the 
obtained F^ t . 

In Fig. [7(a-c) the exact computational results for the lateral Casimir force as a function 
of the phase shift between the sinusoidally corrugations on both surfaces are shown by the 
solid lines at separations a = 124.7, 128.6, and 149.8 nm, respectively (the the first set of 
measurements with corrugation amplitude on the sphere of 13.7 nm). No fitting parameters 
are used in the comparison As can be seen in Fig. [T^a-c), the solid lines are in a very 
good agreement with the experimental data shown as dots. These lines clearly demonstrate 
deviations from the sinusoidal behavior which decreases with the increase of separation. 
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Thus, both the experimental data and the exact theory confirm the prediction of Ref. 
made using the PFA approach that the lateral Casimir force is asymmetric. 

In Fig. [8] the exact computational results for the max \ F^ t \ are plotted as a band between 
the solid lines versus the separation a between the mean levels of corrugations. The width 
of the band takes into account the computational errors equal to 1.5% and the correction to 
Eq. ( 166]) due to surface roughness. The lower solid line is obtained as the computed max \ F^ t \ 
minus 0.015 max \F^ t \. The upper solid line represents (1 + 0.015 + ^ cor r) max \F^ t \, where 
the correction factor 77 corr for the first set of measurement was found in Sec. V. As can be 
seen in Fig. [HI the exact theory is in a very good agreement with the measurement data 
although it deviates significantly from the theoretical band between the two dashed lines 
computed using the PFA approach. Again no fitting parameters are used in the theory. This 
allows one to conclude that with the grating period of A = 574.7 nm the diffraction-type 
effects on the lateral Casimir force have been reliably demonstrated both experimentally 
and theoretically. 

The results of similar computations using Eq. ( 166]) are presented by the solid lines in 
Fig. |9]^a-c) where the lateral Casimir force is plotted as a function of the phase shift for the 
second set of our mesurements at a = 134, 156.5, and 179 nm, respectively. As before no 
fitting parameters are used. In spite of the fact that the experimental data of this set of 
measurements with deeper corrugations on the sphere are more noisy, the exact theory is 
clearly consistent with data and confirms the asymmetric (nonsinusoidal) character of the 
lateral Casimir force. In Fig. [10] the exact computational results for the max|i^ t | in the 
second set of our measurements versus separation are presented as the band between the 
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two solid lines. The width of this band takes into account both the computational errors 
and the surface roughness. It was found using the same procedure, as for the first set of 
measurements, with the replacement of ry corr by fj corr (see Sec. V). As in previous comparisons, 
no fitting parameters are used. It can be seen in Fig. [TD]that in spite of the larger noise in the 
second set of measurements, the data are consistent with the exact theory. The deviations 
of both the exact theory and the measurement data from the theoretical prediction of the 
PFA approach are even more pronounced than in the first set of measurements. Thus, the 
second set of our measurements confirms the observation of the diffraction effects in the 
lateral Casimir force. 



VII. CONCLUSIONS AND DISCUSSION 

In the foregoing we have presented the results of an experimental and theoretical inves- 
tigation of the lateral Casimir force which arises between Au-coated aligned sinusoidally 
corrugated surfaces of a sphere and a plate spaced in close proximity to each other. The 
most distinctive feature of this experiment in comparison with the one performed previously 
371 . |38| is the use of the grating with much smaller corrugation period which permitted to 
enter the region of parameters outside the applicability of the PFA and observe the nontriv- 
ial diffraction-type effects. Another distinctive feature of this experiment is the use of much 
deeper corrugations on both the plate and the sphere. This allowed us to observe an asym- 
metry of the lateral Casimir force (i.e., the deviation of its profile from perfect sine) which 



was predicted [38[ in the framework of the PFA approach, but not observed up to date due 
to the use of small corrugation amplitudes. In this paper we have presented details of the 
experimental setup, the original procedure allowing to imprint relatively deep corrugations 
using a metallized template grating on a sphere of about 100 fim radius, and the process of 
electrostatic calibrations giving the possibility to measure small forces and at short absolute 
separations with high precision. The dependence of the lateral Casimir force on the phase 
shift between the corrugations on the sphere and the plate was measured at many different 
separation distances in the two sets of measurements with different amplitudes of corruga- 
tions on the sphere. The maximum value of the lateral force as a function of separation was 
also investigated in the two sets of measurements. Both the random and systematic errors 
of the measured forces and separations were found and the total experimental errors were 
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calculated at a 95% confidence level. 

The experimental data were compared with an approximate theory using the PEA and 
with the exact theory based on the scattering approach. Note that exact calculation methods 
allowing the evaluation of the Casimir force between bodies of arbitrary shape and made of 



any real material have been actively investigated since 2006 35|, |36|, l6lH7lj. However, until 



the present paper and the rapid communication [41( there was no detailed comparison of 
exact theoretical results beyond the Lifshitz theory taking into account such experimental 
conditions as real material properties of the test bodies, surface roughness and nonzero 
temperature with measurement data of any specific experiment. Using the exact theory, the 
computed dependences of the lateral Casimir force on the phase shift and of the maximum 
magnitude of the lateral force as a function of separation were found to be in agreement with 
the experimental data. Regarding the PEA approach, the computational results were found 
to be excluded by the data. This provides the quantitative confirmation for the fact that 
this experiment was performed outside the applicability region of the PFA and marks the 
beginning of investigations of diffraction effects in the phenomenon of the lateral Casimir 
force. 

As was already noted in the Introduction, the lateral Casimir force might find applica- 
tions in nanotechnology. Specifically, it was proposed to use this effect for the frictionless 
transmission of lateral motion by means of a nanoscale rack and pinion without intermesh- 



ing cogs [28 



, |29[ or in a ratchet with asymmetric corrugations driven by the Casimir force 



30f | . The above experimental and theoretical results bring us closer to the realization of 



such kinds of micromachines. In future it is planned to perform similar experiments with 
corrugations of more complicated shape which provides an opportunity to modulate the 
respective lateral Casimir force. 
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Figures 
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FIG. 1: Schematic of the experimental setup (see text for further details). Insertion shows the 
obstructed side of the sphere with imprinted corrugations. 
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FIG. 2: (Color online) (a) An AFM scan of the grating surface showing the sinusoidal corrugations 
covered with stochastic roughness, (b) A typical section of the grating surface along a y = const 
plane. The solid line shows a sine function obtained from the fit. 
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FIG. 3: (Color online) (a) An AFM scan of the surface of the sphere used in the first set of 
measurements showing the imprinted sinusoidal corrugations covered with stochastic roughness, 
(b) A typical section of the grating surface along a y = const plane. The solid line shows a sine 
function obtained from the fit. 
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FIG. 4: (Color online) The imprinted corrugations on the sphere used in the first set of measure- 
ments. The lighter area shows higher points and hence demonstrates the sphericity of the imprinted 
surface. 
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FIG. 5: (Color online) The imprinted corrugations on the sphere used in the second set of mea- 
surements. The lighter area shows higher points and hence demonstrates the sphericity of the 
imprinted surface. 
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FIG. 6: The deflection signal due to the lateral electric force versus the phase shift (a) before and 
(b) after the correction to the tilt between the grating and the x-axis was introduced. 
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FIG. 7: (Color online) The experiment (dots) and the exact theory (solid lines) for the lateral 
Casimir force versus the lateral displacement normalized for the corrugation period at the separa- 
tions (a) 124.7 nm, (b) 128.6 nm, and (c) 149.8 nm in the first set of measurements with corrugation 
amplitude on the sphere of 13.7 nm. No fitting parameters are used. 
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FIG. 8: The experimental data (crosses) and theoretical values computed using the exact theory 
(the band between the solid lines) and using the proximity force approximation (the band between 
the dashed lines) for the maximum magnitudes of the lateral Casimir force versus separation in 
the first set of measurements. No fitting parameters are used. 
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FIG. 9: (Color online) The experiment (dots) and the exact theory (solid lines) for the lateral 
Casimir force versus the lateral displacement normalized for the corrugation period at the separa- 
tions (a) 134 nm, (b) 156.9 nm, and (c) 179 nm in the second set of measurements with corrugation 
amplitude on the sphere of 25.5 nm. No fitting parameters are used. 
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FIG. 10: The experimental data (crosses) and theoretical values computed using the exact theory 
(the band between the solid lines) and using the proximity force approximation (the band between 
the dashed lines) for the maximum magnitudes of the lateral Casimir force versus separation in 
the second set of measurements. No fitting parameters are used. 
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FIG. 11: (Color online) Two longitudinal gratings of different shape and same period A. 
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FIG. 12: (Color online) The fictitious longitudinal grating for which one evaluates the reflection 

(2) 

matrix R^ t (see text for further discussion). 
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FIG. 13: (Color online) The top longitudinal grating in Fig.[TT]for which one evaluates the reflection 

(2) 

matrix R^' - The normal and lateral shifts with respect to the fictitious grating shown in Fig. [T2] 
are denoted by L and xq, respectively. 
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TABLE I: The mean values (column 2), the variances of the mean (column 3), the systematic 
errors (column 4), and the total experimental errors at a 95% confidence level (column 5) of the 
measured maximum magnitudes of the lateral Casimir force / = max \F^ t \ at different separations 
(column 1) for the first set of measurements. 
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TABLE II: The mean values (column 2), the variances of the mean (column 3), the systematic 
errors (column 4), and the total experimental errors at a 95% confidence level (column 5) of the 
measured maximum magnitudes of the lateral Casimir force / = max \F^ t \ at different separations 
(column 1) for the second set of measurements. 
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